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1. Introduction 

As it is well-known, pure non-abelian gauge theory at T = and space-time dimension d < 4 
possesses the extraordinary property of being in a single phase for all values of the bare gauge 
coupling, < /3 < oo. The theory exhibits passage from the short distance perturbative (ordered) 
regime to a long distance confining (disordered) regime without encountering a phase transition. 
Such multi-scale problems are notoriously difficult to treat from first principles. Another well- 
known, and in fact somewhat similar example is the passage from laminar to turbulent flow in 
hydrodynamics. 

In principle such problems can be treated by some systematic block-spinning or decimation 
procedure. Exact block-spinning schemes in gauge theories, however, so far appear intractable, 
both analytically and numerically. There are, however, approximate but computable decimation 
procedures that can provide bounds on judicially chosen quantities. Such bounds can then be used 
to rigorously constrain the flow behavior of the exact theory [JTJ| - [Q]. The basic procedure followed 
in this program is as follows: 

(i) Devise decimation schemes such that at each decimation step one can obtain both an upper 
and a lower bound for the partition function. Furthermore, the decimations should be explicitly 
computable to any accuracy. 

(ii) At each step interpolate between the upper and lower bound. Then the exact value of 
partition partition function (free energy) is obtained at some particular value of the interpolating 
parameter. This then gives a representation of the partition function on the decimated lattice [[I]]. 

(iii) The same procedure is then applied to the partition function in the presence of a 'twist' 
(external center flux). 

(iv) Combining (ii) and (iii) allows consideration of appropriate long distance order parame- 
ters: the vortex free energy, Wilson loop, and Wilson lines. 

We will see that the flow into the IR regime of the exact theory can be extracted without 
knowing the precise right numerical values of the extrapolation parameters that give the exact 
representations at each decimation step. For SU (2), the only case explicitly considered here, the 
result is flow into the strong coupling confinement regime for any initial j8 . Furthermore, if these 
right numerical values of the extrapolation parameters can be explicitly obtained, one can compute 
appropriate quantities directly on coarser lattices. 

2. Decimation and interpolation 

Starting with some reflection positive plaquette action at lattice spacing a, e.g. the Wilson 
action A p (U) = ^ RetrL^, consider the character expansion of the action exponential: 



F(U,a)=e A ^=^d j F0,a) Xj (U) 



(2.1) 



j 



For SU(2), one has j = 0, \ , 1 , | , . . . 



, dj = (2j + 1). It is convenient to work in terms of normalized 



coefficients; 



F(U,a)=F l + Y,Cj(p)djZj(U) = F f(U,a). 



(2.2) 
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The partition function on lattice A is then 

Z A (J8) = / dU A l\F p (U,a) =F(j A| f dU A Hf p (U,a) . (2.3) 

p p 

We employ decimation schemes of the 'potential moving' type. The lattice is split into hy- 
percubes of side length Xa, and plaquettes are moved from the interior of each hypercube to its 
boundary, where plaquette interactions are appropriately adjusted to compensate for the move. The 
process can be summarized as an explicit decimation rule for each step in successive decimations 

a^Xa^X 2 a^ > X n a <=► A -> A® -» >A (,,) , 

which gives explicit expressions for the computation of the Fourier coefficients at the (m+ l)-th 
step given those of the m-th step: 

F (m) —> F (m + l) = F (m,£,r,X,{ci(m)}) 

Cj(m) -> cj(m + l) = Cj(m,C,r,X,{ci(m)}) 

The rule for the decimation scheme we adopt, which need not be explicitly given here, involves 
parameters £, r which control the amount by which the interactions of the remaining plaquettes 
after a decimation step are renormalized to compensate for the ones that were removed. 

The resulting PF after n decimation steps is then: 

Za(M = fl FoH^"" 1 f dU Ain] Hf p (U,n) . (2.4) 

It is important to note that, after each decimation step, the resulting action retains the original 
one-plaquette form but will, in general, contain all representations: 

f p (U,n) = \l + Y,Cj(n)djZj(U)] =expA p (n) (2.5) 

Mo 

with 

M") = LPj(P>")Zj(U). (2.6) 

j 

Also, both positive and negative effective couplings jSy will, in general, occur. 

The following two basic statements can now be proved. Going from the n — 1 step to the n 
decimation step one has: 

(I) with f = X d - 2 , r = 1 - e , < e < 1 : 

Za(P,»-i)<Za(M; (2-7) 

(II) with C = 1, r=l+£, 0<£: 

Z A (p,n)<Z A (p,n-\) . (2.8) 
In fact, in both (I), (II) one has strict inequality at any finite j3. 
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Translation invariance, convexity of the free energy, and reflection positivity (positivity of 
Fourier coefficients) underlie (I), (II). 

We next interpolate between the upper and lower bounds. Introducing a parameter a, (0 < a), 
define interpolating coefficients dj(n,a) and F (n,a) such that 



Cj (n,t; = X d - 2 ,r= 1-e) : a = 1 
Cj(n, £ = l,r = 1 + e) :a = 



cj(n,a) 
and similarly 

Fo(n,a) = \ F ^ C = ^V = l-e):« = l 
uv ; |F (n,C = l,r=l + e) : a = 

for < £ < 1. A variety of convenient explicit examples of such interpolating cj(n,a), Fo(n,Gt) 
may be given. There is, of course, nothing unique about any one such smooth interpolation. Con- 
sider more generally a family of smooth interpolations parametrized by a parameter t in some 
interval {t\,t%): 

cj(n,a,t) , F (n,a,t) 
Define the corresponding interpolating partition function Z A (fi,n,a,t) by the replacement 

Fo(rc) — *Fo(n,a,t) , Cj(n) — > dj(n,a,t) 
Then, for each fixed t, we have 

Z A (p,n,0,t) <Z A (fi,n- 1)< Z A (p,n,l,t) . (2.9) 
But then, by continuity, there exist a value 

0<a = a K n \t) < 1 

such that 

2 A (fi,n,a£\t),t)=Z A (P,n-l) . (2.10) 

(n) 

In other words, a A (j3,t) is the level surface (implicit function) solution to 

Z A (P,n,a,t)=Z A (P,n-l) . (2.11) 

There is thus 'parametrization invariance' in dnc| ): shifts in t are compensated by changes in 
a^\f},t). This procedure of fixing a^\j5,t) is iterated in successive decimations (n = 1,2,. . .). 
So starting at the original spacing a, after the 1st decimation step the procedure gives an exact 
representation on the coarser lattice of the original partition function on A in the form: 



Z A (J8) =Fj A| JdU A Y[f P (U,a) 



p 

(i) 



z A G3,<(0,0 

F^FoiaW^tp^ J dU Ail) Uf p (U,a£\t),t) . (2.12) 



Since the resulting action in ( |2.12| ) is again of the same type ( |2.5| ), iteration of the procedure results 
in a representation analogous to (2.12) on successively coarser lattices A^. 

Note that, by construction, the coefficients cj, Fo evaluated at a A (t) occuning in (2.12) are 
always sandwiched between the explicitly computable upper bound (a = 1) and lower bound (a = 
0) coefficients. 
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3. Twisted partition function - Vortex free-energy 



The twisted partition function Z A (/}) is defined as the partition function with 'twisted' action 
on each plaquette of a coclosed set Y of plaquettes winding around the periodic lattice in (d — 2) 
directions perpendicular to a given 2-plane. Thus on the dual lattice Y forms a closed 2-dim 
surface in d = 4, and a closed loop in d = 3. The 'twist' is a shift by a non-trivial element z G 
Z(N) of the group center, and amounts to a discontinuous (singular) gauge transformation on the 
configurations in the partition function with multivaluedness in Z(N), i.e. the introduction of a 
lh(SU(N)/Z(N)) = Z(N) vortex. X = -1 for 51/(2). 



The vortex free-energy order parameter F A ' is then defined as 



exp 



,(-) 



^A 

Z A 



(3.1) 



The behavior of the Wilson loop and Wilson line correlations can be related to (3. 1 ) through known 
inequalities. 

There is a slight technical complication in applying the above procedure to the twisted partition 
function. For Z A ' reflection positivity holds only in planes perpendicular to the directions in which 
the set y carrying the twist is winding around the lattice. One way to overcome this is to simply 
replace Z A ' by the quantity Za + Z A ' for which reflection positivity, as it is easily verified, is 
restored in all planes. The above development then can be carried through in this case also. 



Thus, e.g., after one decimation, one has the exact representation analogous to (2.12): 



Z A (P)+Z { - ) (J3)=zJp,a? 1) ( 



t),t)+Z 



(-) 



t)J 



(3.2) 



One may then iterate the procedure in successive decimation steps. 



It should be noted that, as indicated by the notation, the at (t) in (3.2) are fixed indepen- 



dently of and hence a priori can be distinct from the a^\t) occurring in the representation ( |2.12 ) 



for the partition function Za(/3). The same of course is true for each subsequent iteration. 

Consider now using these exact representations for computing the vortex free-energy (|3~l|). 



After one decimation, using (3.2), one has 



Z A + Z 



(-) 
A 



Z A (/3,a+ (1) (0 



t)+z 



j8,a 



+(1), 



t)J 



Z A (a^ l \t),t 



1 + - 



Z H 

^A 



a 



+(i) 



(t),t 



Z A a 



*A 



t),t 



(3.3) 



(3.4) 



By construction, the numerator in ( p.3| ) is invariant under shifts in t; but the two factors in which it 

*A +( 



is split in (^4) are not individually invariant. (They would be, by ( 2.12| ), only if a^ (t) were to 



exactly coincide with a! 1 ' (t) for all t.) 
One now finds that there exist a value t 



t* at which Za( 0f A (t),t ) intersects the level curve 



(/),/ 



i.e. 



Z A [a + A {{ \t*)X 



Z A 
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Thus 

l + ^- = i+ : / +m (3-5) 



$ may then be iterated in successive decimation steps. Now, by construction, the coefficients 
Cj(n, at (t)), Fo(n,a^" \t),t) are sandwiched between the upper bound (a = 1) and lower 
bound (a = 0) coefficients which, for the potential-moving type of decimations employed here, 
are explicitly computable. In particular, the upper bound ones are essentially those of a Migdal- 



Kadanoff decimation. For SU(2) this implies flow of the effective action in (3.5) into the strong 
coupling regime, hence confining behavior for the order parameter, for any initial j8. Decima- 
tions need only be earned out till the strong coupling regime is reached. From that point onward, 
computations can be performed within the convergent strong coupling cluster expansion. 



4. Conclusion - Outlook 



The iterated procedure of decimation followed by interpolation between computable upper and 
lower bounds allows one to constrain and hence follow the qualitative RG flow of the exact theory. 
This is possible without knowledge of the precise numerical values of the extrapolation parameters 
(the Gt^'s) needed to give the exact representation of the physical quantity (partition functions, 
vortex free energy, etc) considered. For SU(2), the case explicitly treated, this suffices to extract 
IR flow into the strong coupling confinement regime for any initial j8 . 

The derivations above only ascertain, within upper and lower bounds, the existence of appro- 
priate values (the a^'s) for the interpolating parameters reproducing a given physical quantity on 
successively coarser lattices. If these exact numerical values of the extrapolation parameters can 
be obtained, one can compute this quantity directly on coarser lattices. An advantage of the above 
development is that the general form of the effective action entering at each step is determined (eq. 
(jO))), and only the numerical coefficients (effective couplings) are to be determined. Algorithms 
for doing this, involving an improved microcanonical demon technique, are under development in 
collaboration with A. Velytsky. 
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